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Abstract 

We study the quantum dynamics of a charged particle in a two-dimensional lattice, subject to 
pH ' constant and homogeneous electric and magnetic fields. We find that different regimes characterize 

these motions, depending on a combination of conditions, corresponding to weak and strong electric 
field intensities, rational or irrational directions of the electric field with respect to the lattice, and 
small or large values of the magnetic (Peierls) phase. 

c3 ' 

1 Introduction and statement of results 

This is the second paper in a row on the quantum dynamics of a charged particle in a two-dimensional 
square lattice, under the influence of an in-plane electric field and a normal to the plane magnetic 
field, both uniform in space and constant in time. In a previous work [T] we considered the case 
where the electric field is aligned with one of the lattice axes. We now extend this investigation to 
arbitrary directions of the electric field, still lying on the lattice plane. New phenomena are found in 
this generalization. 

As discussed in [T], many physical systems, including Bose condensates in optical lattices, can be 
described by the Hamiltonian H of a charged particle in a two-dimensional lattice oriented as the 
reference frame (x, y), in the tight-binding approximation: 
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{m) hm = -y (e- i27ram ^+i,m + e* 2 ™ m ^i, m ) - ^ (Vi, m+ i + 4>Un-i) + ea(FJ + F v m)ip m ,i , (1) 

which is written here for the Landau gauge A = B(—y,0). In the above, tp is the wave-function, the 
integer pair (l,m) labels the lattice site (x,y) = (la, ma), a is the lattice period, J x , J y are the hopping 
matrix elements, F x , F y are the electric field components, e is the charge, and a is the Peierls phase, 
denned as the ratio between the magnetic flux through the unit cell, Ba 2 , and the elementary flux hc/e. 
A second ratio is crucial in this Hamiltonian, that between the Bloch frequencies lo x = eaF x /H and 
ujy = eaFy/h associated with the two components of the electric field vector. For convenience we call 
the ratio /3 = uj x /uiy the (electric) field orientation. These two ratios deserve to be clearly marked for 
further use: 

a := eBa 2 /hc, /3 := w x /w y = F x /F y . (2) 

The content of this paper is a detailed study of the dynamics generated by the above Hamiltonian. 
These dynamics are far from trivial even when one of the fields is absent. In fact, for vanishing magnetic 
field the wave-packet motion features Bloch oscillations, which crucially depend on the rationality of 
the parameter f3. In the opposite case of null electric field, a commensurability condition still appears, 
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this time in terms of the Peierls phase a. Clearly, the problem becomes even subtler when both fields 
are present. To some extent, this problem has been considered earlier in Ref. [H [31 HJ [5]: our goal 
is now to develop a systematic analysis that covers all dynamical regimes. In the paper [1] we have 
developed this analysis in the semiclassical region \a\ <C 1, for the particular case when the electric field 
is aligned with one of crystallographic axes of a square lattice: that is, we examined the case (3 = 0. 
In this work we extend these studies to arbitrary directions of the electric field vector and to arbitrary 
magnetic field intensities. This yields a three parameters space (F, a, (3) (F is the amplitude of the 
electric field) characterized by various dynamical behaviors. We elect not to vary a further parameter, 
the ratio between couplings in the two orthogonal directions of the lattice, J x and J y , keeping them 
equal to J 0. Also, without any loss of generality, we assume |a| < 1/2 and < f3 < 1. 

Our results can be summarized as follows. In sect. 3 we first focus on the case of small a, where we 
use a kind of semiclassical approach in which the Peierls phase a plays the role of an effective Planck 
constant. We derive the one-dimensional classical Hamiltonian 

eaF 

H cl = - J y cos P - J x cos Y + F X P + T y Y , F x , v = (3) 

lira 

as an approximation of the Hamiltonian (flj . Observe that the dynamics of this system is insensitive of 
rationality of the parameter a. We find that both classical and quantum dynamics strongly depends on 
whether the following conditions hold: 

T x < J x and F y < J y . (4) 

First of all, when the above are verified, a stable island gives rise to classical streaming across the lattice, 
in a direction orthogonal to the electric field, with velocity v* given by 

v* = ea 2 F/ha. (5) 

Since semiclassical theory is applicable whenever the Peierls phase a is much smaller than one, we 
therefore predict that when 

eaF/2irJ < a < 1 (6) 

the quantum motion generated by the Hamiltonian (TTJ) allows for wave-packets traveling at speed v* 
for any value of the direction ratio /3, at least for a finite time-span that increases as a tends to zero. 
Conditions © define the small field, semiclassical regime. 

We then turn to a purely quantum analysis. In Sec. 0] we start by examining the case of rational 
orientation /3 and we argue that for non-zero electric field the spectrum of ([T]) is absolutely continuous 
for any value of a. We describe two techniques to compute this spectrum, that can also be implemented 
numerically. Using spectral analysis, when (|6]) holds, we construct the Stark transporting states, already 
defined in the particular case Ref. [1] , that quantize the semiclassical transporting islands of Hamiltonian 
([3]). We show that their evolution is characterized by linear motion with speed v*. Also in Sec.SJ we use 
perturbation theory to compute the leading order term, in the inverse field amplitude 1/F, of the width 
of the energy bands of Hamiltonian (TTJ). Letting f3 — r/q, with r,q co-prime integers, we find that these 
widths scale as (1/ F) q+r ~ l . This suggests that the energy spectrum of this Hamiltonian is pure point 
for irrational direction ratio j3. Note that being based on quantum perturbation theory, these results are 
not restricted to the semiclassical region. 

Next, we turn to the numerical simulation of the time evolution of an initial wave-packet, to confirm 
the theoretical predictions of the spectral analysis. The numerical techniques are briefly described in 
Sec. [5] They permit to attack the most general case, that includes non-rational directions of the electric 
field. First we consider the case F ^> F cr = 2ira.J/ea, where we can apply the perturbation theory of 
Sec. HI and show that an initial wave-packet spreads in the direction orthogonal to the electric field with 
a dispersion a that increases in time as 

a „ (l/F^+^h (7) 
(here, as above, /3 = r/q). This behavior sets in after a transient time, 

T-F^+r-V, (8) 
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during which the packet is localized. Thus for irrational field directions the wave packet is always 
localized in the lattice, which is consistent with discreteness of the spectrum for irrational ft. 

Next, in Sec. [6] we describe the case F < F cr for the semiclassical region \a\ <C 1. We find that the 
motion is either a ballistic spreading (for an incoherent initial wave-packet) or directed transport (for a 
properly devised coherent initial wave-packet), both for rational and irrational directions ft, for all times 
that we could reach in our numerical simulations. This seems to be inconsistent with discreteness of 
the spectrum for irrational ft. To resolve this seemingly contradiction we turn in Sec. [7] to the case of 
arbitrary values of the Peierls phase a. We simulate the quantum evolution for an increasing sequence 
of values of a while keeping the classical parameters T x and T y fixed. For irrational directions ft, at 
fixed time, we observe a sharp suppression of ballistic spreading, occurring when a overcomes a certain 
threshold. Turning then to a time-resolved analysis, we observe saturation of the ballistic spreading, 
for times larger than a threshold that grows abruptly with the inverse of the electric field intensity F. 
Combined with the observation on band widths, this fact strengthen the conjecture that the spectrum 
of the system is pure point, for any irrational direction ft; yet the localization length should depend 
sensitively on the magnitude of the electric field F, like in certain two-dimensional models of quantum 
rotators and of Anderson localization [8j |9j [10] . 

2 Aligned electric field: review of results 

To introduce concepts and notations, we now briefly review the results of our previous work pQ, where 
we considered the case of an electric field F = (0, F) aligned with the y axis of a reference frame defined 
by a two-dimensional square lattice. The magnetic field is directed orthogonally to the lattice plane, and 
the Peierls phase a is small. The system shows two qualitatively different dynamical regimes, depending 
on the magnitude of the electric field: a regime of directed transport for weak fields and a regime of 
ballistic spreading for strong fields, separated by a critical magnitude, F cr = 2iraJ x /ea. Although we 
will keep track of physical constants for completeness in what follows, in numerical experiments we adopt 
adimensional units so that e = a = H = 1 and we set J x = J y = 1. Also, until Sec. [7] or unless otherwise 
noted, the Peierls phase is a — 1/10. 

2.1 The energy spectrum 

The stationary Schrodinger equation derived from Eq. (TTJ) is 

- y (e- l2 ™ m ^+i,m + e i27rQm ^-i,m) - y 1 (VWl + ^i.m-i) + eaFmip l>m = E^ m . (9) 

Imposing a periodicity of period L x = La in the /—direction in Eq. @, where L eventually tends 
to infinity, and using the substitution ipi tm — b m (n), where n is the dimensional quasimomentum, 
k = 2nk/aL, with k = 0, 1, . . . , L — 1, taking values in [0, 2ir/a), one reduces ((9]) to the following ID 
equation for the coefficients b m (n): 

- y (b m +i + b m -i) - J x cos(27ram - an)b m + eaFmb m = Eb m . (10) 

The spectrum of (|I0[) for a fixed k is a modulated Wannier-Stark ladder and the energy bands appear 
by varying the quasimomentum. For example and for the purpose of future comparison, the left panel of 
Fig. [U shows the eigenvalues E v (n) of this equation for F = 0.3. To understand the displayed spectrum 
the semiclassical approach proves to be very useful. 

2.2 Semiclassical theory 

When \a\ is much less than one, the period of the oscillating phase in Eq. (|I0[) is much larger than the 
lattice period, so that we can approximate the discrete function b m by a continuous function B of the 
coordinate y. Furthermore, consider the shift operator T a = exp(ad y ), that acts as T a B(y) = B(y + a). 
Using this in Eq. pH|) leads to 

[— J y cos(iad y ) — J x cos(27ray / a — an) + eFy]B(y) = EB{y). 
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Figure 1: (Color online) The energy spectrum E = E v (k) for F — 0.3 and (r,q) = (0,1) (left) and 
(r, q) = (1, 1) (right). The vertical line in the right panel marks the first Brillouin zone, the size of which 
scales as 1 / yfr 2 + q 2 . 

Next, introducing the operators Y = 2-Kayja and P = —iad y in the above leads us to the effective 
Hamiltonian H qu = — J y cosP — J x cos(Y — an) + J-Y, whose classical counterpart reads 

= -J y cosP- J x cosY + FY, F=^-. (11) 

Since in the quantum description the canonical variables P, Y are operators obeying the commutation 
relation [Y, P] — i2ira, the semiclassical parameter is here the Peierls phase a. If \a\ -C 1, one can 
explain certain features of the spectrum in Fig. 1(a) by simply analyzing the phase portraits of the 
classical system (|11[) . In particular, it is easy to see that the Hamiltonian (|lip can support bounded 
motions if F < F cr : these are seen as the islands around the fixed points in Fig. [21 that draws the 
solutions of Hamilton equations for the Hamiltonian H c i, Eq. (fTTj) . 

By introducing the Bloch frequency u = eaF/h and using the canonical substitution P — > P — uit, 
one obtains from (jlip the time-dependent Hamiltonian 

H c i{t) = -J' y cos{P -ojt)- J' x cosY , J' Xiy = 2naJ Xt y , (12) 

in which the canonical variables lie in the torus, — ir < P,Y < tt. The classical phase space of solutions of 
(fTTI) , pictured in Fig. can also be calculated as the stroboscopic map of the time-periodic Hamiltonian 
(|12l) . Two cases for F < F cr and F > F cr are shown in Fig. [^a,b). Bounded trajectories of pT|) appear 
now as two nonlinear resonances (transporting islands), whose size shrinks to zero when F = F cr . 
Quantizing these islands leads to the transporting modes, defined below. 

Note that in Ref. p] a different gauge, A = B(0, x), was employed. Needless to say, the quantum 
spectrum of the system does not depend on the gauge adopted. In classical mechanics, the latter gauge 
corresponds to the Hamiltonian H c i{t) = —J' x cosP— J' y cos(X + tot), that can be obtained from (fT2")l by 
a canonical transformation that interchanges the role of P as momentum and Y as coordinate. 

2.3 Directed transport regime 

When F is small, the energy bands E(k) overlap forming a rather complicated band pattern. If F < F cr 
this pattern contains a number of parallel straight lines, seen in Fig.[TJa). In the semiclassical approach 
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Figure 2: Phase portrait of the classical system (fT4f for (r,q) — (0,1) (upper row) and (r, q) = (1,1) 
(lower row) at electric field intensities i* 1 = 0.3 (left column) and F = 1 (right column). Twenty different 
trajectories are shown in each case. 



these lines are associated with nonlinear resonances, that transport the particle in the x direction at the 
drift velocity ([5]). The explicit form of the quantum transporting states is 



g(K)e lKal b m ( K )d K , 



(13) 



where b(«) is the eigenvector of (1101) associated with a given straight line in the energy spectrum and 
g(k) is an arbitrary square integrable function of the quasimomentum. Note that in b(«) we jump over 
avoided crossings when following a straight line in E(n): this is the so-called diabatic approximation. 
Thus the vector function b(«) in Eq. fjl3|) . which implicitly assumes the extended Brillouine zone picture, 
satisfies the translational relation 



b m (n + 2n/a) rj & m +i/ a («) , 

which is exact if 1/a is an integer. This relation should be compared with the relation b m (n + 2n /a) = 
b m (n) which holds when we follow the energy bands adiabatically. 

The left row in Fig. [3]shows transporting states for g(n) ~ exp[— C(aK/2n) 2 ] for two different values 
of the parameter C [TT] . Since the slope of the straight lines in the spectrum is equal to the drift velocity, 
the time evolution of these states is a shift of the wave packet in the positive direction at velocity v*. 
Thus the wave-packet width a is constant in time, while the first moment grows as M\{t) — v*t. Here and 
below we identify the wave-packet width with the dispersion, i.e., a — y/ M% — Mf, where M», i = 1,2 
are the first two moments of the position operator. 

It is worth stressing that the above statement refers only to initial conditions given in Eq. (|13[) . For 
a generic initial wave packet one observes an asymmetric ballistic spreading (see Fig. 10 in Ref. PQ), 
where M\(t) w and a(t) — > At. The coefficient A is here defined by the drift velocity: A ~ v* ~ F and 
not by the hopping matrix element J x as in the strong field regime described in the next subsection. 
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Figure 3: Example of transporting states IvE^ml 2 for (r, q) = (0,1) (left column) and (r, q) — (1,1) 
(right column) depicted as grey tone images (black maximum, linear intensity scale). The upper row is 
characterized by the width parameter (7 = 1, while the lower row has C — 0.2. The other parameters 
are a = 0.1 and F = 0.1. 



2.4 Ballistic spreading regime 

Consider the strong field limit, eaF ^> J x ,y Using first order perturbation theory the solution of (JTUJ) 
reads 

E v {k) = eaFv — J x cos(<zk + 2irav) , 
b m\ K ) = <W ± {Jy/eaF)5 v . m ±i . 

It follows from the above equations that the time evolution of a generic localized wave packet is a 
ballistic spreading in the positive and negative x directions, so that the packet width a grows asymptoti- 
cally as a(t) — At. The numerical factor A in this asymptotic expression depends on the particular form 
of the initial wave packet, yet it is superiorly bounded by A = J x /\^2. The maximal value is reached for 
an incoherent wave packet, that we can mimic by assigning random phases to the complex amplitudes 
of an initial gaussian wave packet. Further down in this paper, Fig. [5] display the F-dependence of the 
coefficient A, obtained by direct numerical simulation of the system dynamics, versus F = F y . The 
coefficient A tends to zero in the limit of weak electric fields, while in the opposite limit of strong fields 
it approaches the constant value A = J x /^/2. We will now investigate how this situation changes when 
the electric field is not aligned with one of the axes of the crystal. 



3 General direction of the electric field: Semiclassical theory 

The classical Hamiltonian ([3]) for generic direction of the electric field, 

en F 

H cl = -J y cos P - J x cos Y + F y Y + F X P , F x , y = , (14) 

Zna 

can be guessed a straightforward generalization of the classical Hamiltonian (fTT|) . More rigorously, the 
same expression can be obtained as done in Appendix, starting from the quantum equation (|21[) below. 
Again, it is convenient to use the canonical transformation P — > P — u> y t and Y — > Y + uj x t, which leads 
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to the time-dependent Hamiltonian 

H c i(t) = -J' y cos(P - uj y t) - J' x cos(r + u x t) , J' Xiy = 2naJ x>y , (15) 

Now the phase space of the system ([15]) can be reduced to the torus only when two Bloch frequencies 
uj x and ujy are commensurate, i.e., if /3 :— uj x /uj y — F x /F y — r/q (r,q are co-prime integers). However, 
semi-classical analysis can be applied equally well in the both cases of rational and irrational f3. Let us 
also note that the system (fT"5|) has the global integral of the motion, 

1 = H cl {t) + F y Y + F X P , 

that reflects the fact that the original Hamiltonian Eq. ^ is one-dimensional and hence is trivially 
integrable, for any value of the ratio /?. 

The lower panels in Fig. [2] show the classical phase space for (3 — 1: classical islands appear for small 
field magnitudes. The structure of phase space changes when the field magnitude is increased, similar to 
what happens in the case /3 = (higher panels) . The transporting islands disappear when at least one of 
the two conditions, T x < J y and T y < J y is violated: this yields the condition (QJ. This result indicates 
that the original quantum system must have two qualitatively different regimes, which are separated 
by a critical field magnitude F cr = 2naJ/ea. This conjecture is supported by analysis of the energy 
spectrum of the quantum system, that we consider now, starting from the case of rational orientation j3. 



4 Energy spectrum for rational orientation of the electric field 



We now proceed to the case of arbitrary, yet rational, direction of the vector F = (F x ,F y ). We first 
review two techniques to compute the energy spectrum in the case of rational ratio (3 — F x /F y . This 
analysis facilitates the understanding of the wave-packet dynamics, that will be described in Sec. [5j 
Letting r and q be relatively prime integer numbers, we align the electric field with the (r, q) direction 
in the plane: 

F = ^=(r,q) (16) 

where N = r 2 + q 2 . We choose the gauge 



A = B - 



N 



q(rx + qy) r(rx + qy) 



<r ' (17) 

which reflects the geometry induced by the electric field. Within this gauge, the Hamiltonian becomes 



(Hi/})i, m = (exp -i2ira^-(rl + qm) 



N 



T \ (vl ~\- QTTl) 

i2ira—{rl + qm) ipi >m+1 + h.c.J + eaF j= — Vi,m ■ 



N 



(18) 



where h.c. denotes the terms required to render the Hamiltonian Hermitian. 

The calculation of the energy spectrum can be equally achieved by cither of two different, yet equiv- 
alent methods, which we shall refer to as the method of rotated coordinate frame, introduced in [2], and 
the method of rotated basis. 



4.1 Rotated coordinate frame 

The method of rotated coordinate frame [2] consists of two steps. The first is to choose the previous 
gauge, Eq. (TlTl) . for the magnetic field. The second step is to simplify the Hamiltonian (| by rotating 
coordinates to align the electric field with the vertical axis £ of a new coordinate frame (f?,£): 

qx~ry rx + qy 

V = i ; , K = i 9 ■ (19) 
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Figure 4: Extended lattice for the field direction (r, g) = (f , f ), where it consist of two sublattices. 

In the rotated coordinates, the original lattice sites (al,am), with integer I and to, appear to lie on a 
sublattice immersed into a new square lattice of spacing d, 

d = a/VN , 

whose sites (sd,pd) can be labeled by pair of integer indexes (s,p). Note that this new lattice actually 
consists of N independent sublattices, only one of which coincides with the original lattice, see Fig. 01 
Explicitly, the mapping of the original lattice points into the new, extended lattice is 

(Z, to) — > (ql — rra, rl + qra) := (s,p), 

while the inverse mapping reads 

(s,p) -> j-f(qp~rs,qs-rp) := (l,m). 

Letting now 4>s,p denote the wave-function amplitude at site (sd,pd) in the rotated frame of reference 
(j], £), and using the fact that (rl + qm) = p, we can write the Hamiltonian action as 

(j?0)., p = -y (e-^^^p+r + h.c) ~ i (j 2 ™ r P' N 4> s - r ,p+ q + h.c) + edFp^ p . (20) 

Observe that, coherently with Eq. (fl8l) . the N sub-lattices described above are uncoupled, and that the 
shift s — s- s + N is an invariant transformation. Nevertheless, it is convenient to solve the eigenvalue 
problem for all sublattices simultaneously. We therefore consider the stationary Schrodinger equation 
for the complex amplitudes <fi, (H<f>) sjp = E<fi StP . Following [5] we use the plane wave basis 

4>s, P = > 

where L eventually tends to infinity. In so doing, the quasi-momentum k belongs to the interval 
[0, —yN). We finally arrive at the following equation, where we have put 9 = 2iraN~ 1 : 

_J_x_ ^ e - t e qP + iq d Kbp+r + e ie q (p-r)-iqd Kbp _^__y ^evp-trd*^ + e -iSHp-<i)+ird Kb \ +edFpbp = Ebp 

2 2 (21) 

Equation (fTU|) in Sec. [5Ja) is the particular case (r,q) = (0, 1) of this more general equation. It is also 
possible to show (see Appendix) that in the semiclassical approach this equation yields the classical 
Hamiltonian (|14[) . after applying a canonical transformation. 
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4.2 Rotated basis 



In this method we do not introduce sub-lattices but diagonalize the Hamiltonian by changing the basis. 
One begins with the tight-binding Hamiltonian in the gauge (fl7|) as before, that is, Eq. (|18p . One 
then imposes periodicity on the lattice space in the following way: let K be an even integer, and let us 
identify lattice points (I, m) and (£', ml) if and only if I — V is an integer multiple of Kq and m — w! is 
an integer multiple of Kr. The Hamiltonian (TT51) typically does not satisfy this periodicity. Nonetheless 
observe that by letting K tend to infinity the period can be made arbitrarily large, following the standard 
approach. At finite K the lattice is then composed of K 2 qr distinct points. 

Next, we define new lattice period 

d=aVN . 



and an adimensional quasi-momentum k — dk, that can take the discrete values 

k = 2wj^, j = 0,...,K-l, (22) 
and a new set of basis functions associated with k: 

-ikn 



I^K> = E ^IP + m^-rn). (23) 

— -K/2 1 1 * 



This definition is consistent with the imposed periodicity on the lattice, as can be seen by letting 
n — » n + K and observing that the r.h.s. does not change. Also observe that if k ^ k' 

<^>K> = o, 

for all choices of p,p' ', /i, p,' . Moreover, keeping k fixed, different functions \4^],) and \<p^\) overlap if 
and only if (p,fi) and (p',/J,') belong to the same discrete transverse line of direction (q,—r): i.e. there 
exists and integer j such that (p', //) = (p, fi) + j(q, — r). In this case, 

The fact that different functions may overlap is clearly a consequence of the overdetermination of 
the set Irf^l), which is made of K 3 rq elements: the integer variables p and /i can take all values from 
to Kq — 1, and from to Kr — 1, respectively, while k can take the K discrete values in Eq. (|2"2"1) . To 
the contrary, as observed above, the periodic lattice requires a basis of K 2 rq elements only. Since for 
each pair (p,fJ>), the function \4>^k) ^ s a li near combination of K lattice sites \l,m), it is easy to see that 
this set can be equally spanned by keeping (p, (i) fixed and letting k vary over the K discrete values in 
Eq. ([2"2")l . It is equally easy to find that the full periodic lattice can then be spanned by choosing (p, fi, k) 
in the discrete set I: 

l = {0,...,Kq-l} x {0,...,r-l} x {0,2tt-^ . . . ,2tt^^}. 

Remark that in the set I the variable /i takes on a restricted set of values, of cardinality r. Equivalently, 
we could have chosen a set organized in strips parallel to the m axis, but since r < q, this choice would 
have been less convenient. To sum up, we have introduced the basis for the periodic lattice that, because 
of the above computations, is composed of orthonormal functions. 

Let us now compute the matrix elements of H over this basis. The fundamental point is that H 
does not couple functions with different fc. In fact, when computing the non-diagonal couplings, the 
product rl + qm appears, that is constant on the transverse lines defined above. The sum over n of the 
phases stemming from Eq. (f2"3"| is then zero, unless fc = k 1 . The same result also holds for the diagonal 
term (F x l + F y m), since (F x ,F y ) is proportional to (r, q), see Eq. (TT6l) . Observe that this is true at any 
finite value of K, independently of the fact that the Hamiltonian H is, or is not, periodic. Therefore, 
the spectral problem of H over the periodic lattice decomposes into the K fibers obtained by letting 
k = 27r j^, with j fixed. At this point one is left with the computation of the matrix elements within a 
fiber. 
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Suppose now that 1 < r < q, we shall treat the remaining cases 1 = r < q. separately. Then, by 
explicit calculation we find that the non-zero matrix elements of H are the following. The "horizontal" 
coupling leads to the matrix element 

(4>Hx,k\H h 4i) = -yexp (W*|(rp + g/ i)) , 
if < p < Kq - 2, while for p = Kq - 1 

<<2l^g-l, fe > = - y exp (»27ra£(r(tfg - 1) + • 
The "vertical" coupling gives 



cxp ^— Vina— (rp + qp,)j , 



if < (i < r — 2, while for fj, = r — 1 



6 (0) I 



1} > = -y ex P (-i27ra£(rp + g(r - 1))) e ifc , 



where p + g at l.h.s. is to be understood modulus Kq. The above are half of the required formulae: we 
must also add the Hermitian conjugate matrix elements. Finally, the diagonal coupling is 

<^|^JS> = ^(rp + gM). 



N 

Observe finally that this lattice problem can be seen as a combination of r coupled, one dimensional, 
periodic lattice problems of size Kq. This takes a particularly simple form in the case (r, q) = (1, q). In 
fact, we end up with the single, one-dimensional lattice of size Kq, on which the fiber Hamiltonian has 
matrix elements 

= "f Wfi + h.c.) - h ( e -^V fe V, p+ , + h.c.) + ^V,p- (24) 



We are therefore equipped with the analytical tools to compute the energy spectrum of the system. 

The right panel in Fig. Q] displays example of the spectrum for (r, q) = (1, 1), which was calculated 
by either of two methods. Note that spectrum is periodic on the quasimomcntum with the period 
2n/d = (2n/a)/\/N. Also note that both spectra in Fig.[T]show the characteristic pattern with straight 
lines: moreover, the slope of lines is the same in both panels. As it will be shown later on, this slope 
defines the velocity of quantum particle in the transporting regime. 



4.3 Continuity properties of the spectrum 

In the preceding paper [l a , both theoretically and in the numerical examples, we have mainly considered 
rational and small values of a, a fact that assured both periodicity of the reduced one-dimensional 
Hamiltonian in (ITU)) and validity of the semiclassical analysis. In this section we discuss the spectral 
properties of the Hamiltonian (UJ) for the case of a general value of a. In this section we consider only 
rational values of the field direction /?. 

As remarked in the introduction, it is instructive to think of the limiting case of null magnetic and 
electric field: the spectrum is absolutely continuous, composed of a single band. Keeping the electric 
field null, while turning on the magnetic field, the one-dimensional Hamiltonian (|10[) . [or equally well 
the Hamiltonian (I2T1) -being the electric field null, the two are equivalent], becomes the celebrated 
Harper Hamiltonian, whose spectrum has bands for a rational: therefore, the spectrum of the original 
two-dimensional Hamiltonians (j9j) and (fT8|) . is also absolutely continuous. For irrational values of a, 
the spectrum of (fTOj) . ([2~T j) is absolutely continuous (for J y > J x ), pure point (for J y < J x ), or singular 
continuous, in the critical case J x — J y . In the first case, obviously, the spectrum of (|5j). (fT5)l is absolutely 
continuous. Interestingly, the same is true also in the second case, J y < J x : by varying the quasi- 
momcntum associated with the invariant direction x the eigenvalues of (|10[) and (1211) . that exist since 
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the spectrum is point, move and draw the energy bands of (|9|) and (|18[) . Finally, for J x = J yi since 
in this case the spectrum of (fTU|) . (|2U is the same for almost all values of k, due to ergodicity of the 
Hamiltonian, the spectrum of (|9j). (118l) is also singular continuous. This is the complete description in 
the case of null electric field. 

The previous results are classical in historical sense. They are useful to understand the non-zero 
electric field case. Since j3 is rational, we can still use the Hamiltonians (fT8 ]) . (|21[) . In this case the 
spectrum of the one-dimensional Hamiltonians (|21[) is pure-point, for any value of a - whether rational 
or irrational, small or large. Therefore, we typically expect the spectrum of (TT5)) to be absolutely 
continuous, drawn by the eigenvalues E v {n). Of course, this is heuristic discussion is not a substitute for 
a rigorous proof, but it can give us a hint into the physical phenomenon. We therefore conjecture that 
for any value of a, and non zero electric field with a rational direction ratio, the spectrum is absolutely 
continuous. Quantum dynamics under these circumstances should reflect the spectral properties and be 
ballistic. We shall prove numerically that this is indeed the case. 

4.4 Strong field limit and perturbation theory 

In the strong held limit we can find analytic approximations of the spectrum by perturbation theory 
based on Eq. (f2"Tj) . These results are independent of the rationality of a and of its amplitude. The case 
j3 = has been described in Sec. 12.41 The next simplest case is (3 — 1, i.e. (r, q) = (1, 1). In this case, it 
is convenient to rewrite Eq. (|21|) in the form 

-[V(p;K)b p+1 + V*(p-l;K)b p - 1 ] + edFpb p = E{n)b p , 

where V(p; re) = ( J x e~™ ap e idK + J y e nav e^ ldK ) /2. Similar to the case (r, g) = (0, 1) analyzed in Sec.jSH 
the unperturbed spectrum of the system consists of flat bands separated by the Stark energy, i.e., 
E®(k) = edFv. However, unlike the case (r,q) = (0,1) the first order correction to this spectrum 
vanishes. The second order correction is given by 

AE U = |V(i/ ~ 1;K)|2 - |y( ^ )|2 = ^-[cos(27ra(i/ - 1) - 2dn) - cos(27rai/ - 2dn)\. 
edF edF 2edF 

This equation proves that the band widths decrease as \/F when F increases. This is a special case of 
a general perturbation theory result. To treat all other cases (r,q), notice that according to Eq. (|2"Tj) , 
the Hamiltonian operator can be written as H — Hq + V , with 

H =J2edFp\p){p\, (25) 
p 

V = ^$(p, g ,K)|p)(p + r| + $*(p-r,q, K )|p)(p-r|+e(p,r, K )|p)(p + g| + e*(p-g,r, K )|p)(p-g|, (26) 
p 

where we have put 

*(p, q, «;) = -y e -^qp e iqdn^ ^ K ) = jOrv e ~irdK _ (37) 

Two properties of this Hamiltonian are immediately noticed. Firstly, as noted above, the unperturbed 
spectrum is equally spaced. Secondly, the perturbation V only couples unperturbed states (here, Kro- 
necker deltas at site p) with quantum numbers differing by either q of r. For short, denote by V p y the 
matrix elements of the perturbation V in the unperturbed basis. 

When q > r > 1 it is then immediate that the diagonal terms of the perturbation V vanish. Proceed 
next to consider second order perturbation theory: it is composed of terms of the kind \V p , p +j\ 2 / (E® — 
Ep + j), where E® is the unperturbed spectrum. In our case, j can only take the values ±g, ±r. Yet, 
for any of these j, the energy difference [E® — E® + j) = —edFj is the opposite of that of the term — j 
(the unperturbed levels are equidistant), so that the two contributions cancel exactly, and second order 
perturbation theory yields a null result. 

The last observation has also a bearing on the general, n-th term in the Rayleigh-Schrodinger pertur- 
bation series for n > 3, that contains the second order term as a factor, among other terms. Typically, 



11 



the n — th term is composed of sums of products, call them JJ, of n matrix elements of the perturbation 
V. The first term in such summation is always of the kind 

n_ Yliil YiLi2l ' ' ' YhirliE. o$\ 
~~ (E° -E°)(E a -E°)---(E° -E° )' 1 1 

Other terms in the summation giving the n-th term contain shorter "chains" of products of matrix 
elements, of length at most n — 1, appropriately multiplied among themselves to give order n, of course 
divided by the related denominators. In some of these, the second order term appears as a factor. 

Notwithstanding this complexity, the particular form of the Hamiltonian problem (|25|27|) permits us 
to derive a simple result. In fact, the numerator in Eq. (|2"5)) is null unless the "path" p j\ — > — > 
jn-i —> P is composed of "allowed" jumps of size ±q, ±r. It is easy to realize that, in order for this path 
to comprise the least number of jumps, it must contain a positive number, s, of steps of length q and 
a negative number, t, of steps of length r (or the same with opposite signs), so that sq + tr — 0. Now, 
since q and r are relatively prime, the minimal solutions are s = r, i = —q or s = —r, t = q. This implies 
that the minimal number of steps must be n = r + q. Therefore, when q > r > 1 the first non-zero term 
in the Rayleigh-Schrodinger perturbation series for the eigenvalues is of order n = p + q and is composed 
of 2(p+q) addenda of the form (|2"5)l . All other terms that appear formally in the analytical expression 
for the n-th term contain shorter "chains" (and/or the second order term) and are therefore null. 

The complete expression of the leading perturbation term can be explicitly computed, using Eq. 
(1271). Yet, notice that all terms of the form (|28|) contain the common factor 



\ , 77j _ ( JxYi J V ) .»,); 

and therefore we can write 



E„( K )~E° + A r , q (F)P n (K), (30) 

where P n (re) is a trigonometric function of the quasi-momentum n. A consequence of this result is that 
the rate of ballistic spreading discussed in Sec. 12.41 is strongly suppressed for large electric field as soon 
as the vector F does not point to a strongly rational direction f3. Later on, in Sec. 16.21 we confirm this 
expectation, as well as the quantitative dynamical estimates that follow from Eqs. (|29I30|) . 



5 Wave packet dynamics: numerical techniques 

To simulate numerically the quantum evolution of the system we have adopted two different approaches. 
In the first, we have used the time-dependent gauge A(t) = Ao + c(F x t, F y t), for which the electric field 
appears as a periodic driving of the system. For the static vector potential Ao, which is responsible for 
the magnetic field, we used the Landau gauge A = £>(0, x). This leads to the following Schrodinger 
equation with explicit time-dependence 

ini>i,m = -y (e-^Vm,m + h.c.) - ^ (e-^^+^^+i + h.c.) , 

that we have solved employing the standard Runge-Kutta techniques implemented in Matlab. 

In the second approach, implemented in Fortran, we have chosen the time independent gauge A = 
B(—y,0), so that the Schrodinger equation reads 

ihi>l,m = - y {e' l27Tam ^i+i, m + e l2 ™"Vz-i,m) - y (Vwi + V^m-i) + ea{FJ + F y m)^ m . 

The r.h.s. of equation is in the form of the action of a time-independent operator on the two-dimensional 
lattice vector ip. We have computed the exponential of this operator using the repeated Chebyshev 
expansion 12 , combined with a numerical truncation of the infinite lattice to a strip along the line of 
equation F x l + F y m = 0, that corresponds to the direction of spreading of the wave-packet. This is 
effected numerically by introducing slanted integer coordinates I' — I and m! = m + mt(-^N). Observe 
that the directions of increase of I' and m! are no more orthogonal. Yet, this provides a convenient 
re-labelling of lattice sites to enforce the strip truncation mentioned above. We always checked that the 
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wave-packet projection at the boundaries of large |/'| and \m'\ is smaller than a very low threshold, at 
any time during the evolution. 

The initial conditions, obviously common to both approaches, consisted of two-dimensional Gaussian 
wave-packets, centered at the origin of the lattice, of adjustable widths: tfji, m ~ exp(— C x l 2 — C y m 2 ). 
To simulate an incoherent wave packet, all components i/'i.m have been multiplied by statistically inde- 
pendent random phases e l1?,>m , and an average of quantum amplitudes over different realizations of the 
initial packet has been performed. 

Finally, to simulate the dynamics of transporting states we have first constructed them by using the 
appropriate gauge for the magnetic field and then we have applied a unitary transformation to translate 
them into the fixed Landau gauge. For example, for the transporting states (|13|) . which were constructed 
using the gauge Ao = B(—y, 0), the unitary transformation reads ipi,m(t = 0) = exp(— i2iralm)'fyi tm . 

6 Wave packet dynamics: Semiclassical region 

Let us start our analysis of the system dynamics from the semiclassical region |a| <C 1. It is instructive 
to have a pictorial look at four significant cases: we select two values of the field intensity, the first weak, 
F = 0.2 and the second strong, F = 0.5. We combine these with two values of the orientation of the 
electric field: a rational value, /3 = 2/3 and a strongly irrational one, f3 = (vo — l)/2. 

Since the wave-packet spreads in the direction of the 77 coordinate of the rotated frame of reference 
in Eq. (|19j) . we find convenient to adopt this coordinate system. Moreover, since we observe that the 
wave-packet is localized in the £ direction, we compute the wave-packet projected amplitude \ip(r],t)\ 2 = 



In Fig. [5] we display \ip(j],t)\ 2 after averaging this quantity over a number of different realizations 
of the random phases in the initial gaussian wave-packet. In the first case, Fig. [5ja), the electric field 
amplitude is F — 0.2, We observe for both orientation of the field sub-packets moving at constant speed, 
in the positive and negative 77 directions. The situation changes radically in Fig. [5jb), drawn on the 
same scale, now for F = 0.5. Here, we find sub-packets travelling in both directions only in the case of 
rational orientation of the field, while in the irrational case we observe localization of the motion also in 
the 77 direction. 

To the contrary, the projected distributions on the £ direction (not reported here), initialized to a 
Gaussian, settle to a shape that is still approximately Gaussian, whose width naturally depends on the 
amplitude of the electric field: a clear sign of Stark localization in the direction of the electric field, that 
does not depend on the orientation of this latter. Let us now consider separately the weak and strong 
field cases. 

6.1 Weak field limit and transporting states 

Suppose now that conditions © hold. Then, according to Sec. [31 F < F cr and the system has transport- 
ing islands. The quantum- mechanical signature of these islands are straight lines in the energy spectrum, 
that are clearly observed in Fig.[TJ The slopes of these lines coincide with the drift velocity v* of Eq. ([5]). 

Adapting Eq. (fT3|) of Sec. 12.21 to the present case of rational field direction, we can construct a 
family of localized wave packets, which move at the drift velocity in the direction orthogonal to the field, 
without changing their shape. To do this, we first write these states on the extended lattice that, we 
recall, consists of N sublattices: 



In the above, <?(«) is an envelope function, that we choose of the form g(n) ~ exp[— C(dn/2ir) 2 ]. Then we 
select from this (generally very large) array only the complex amplitudes ^!^ m which sit on the original 
lattice. We plot these states in Fig. [3] the upper row of this figure displays the transporting states ^i. m 
for (r,q) — (0, 1) and (r, q) = (1, 1), where we choose C = 1. For this value of C the states are equally 
localized in both directions, parallel and orthogonal to the vector F. If C is decreased, the states become 
more extended in the parallel direction and more localized in the orthogonal direction, see the lower row 
in the figure. In the opposite case, i.e. when C is larger than one, the situation is obviously reversed. In 
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Figure 5: (Color online) Wave-packet amplitudes \ip(r),t)\ versus r\ and time, averaged over 12 realiza- 
tions of the initial incoherent wave-function, for J x — J y = 1, a — 1/10 and two electric field orientations: 
(3 = 2/3 (red curve) and j5 = (y/5 — l)/2 (green curve). Electric field amplitude is F — 0.2 (top) and 
F = 0.5 (bottom). 
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Figure 6: (Color online). The center of gravity of the specially constructed wave packet x(t) (blue) and 
y(t) (red) for eaF/J = 0.1 and (3 — 1/3 (solid lines) and (3 — (\/5 — l)/4 (dashed lines) versus time, 
measured in units of the tunneling period period Tj — 2ttJ jfi which, in the units used, is equal to 2w. 
The lower panels show the wave packets at the end of numerical simulations, t — 607r, depicted as grey 
tone images (black maximum, linear intensity scale), for these two cases. 



the limit C —¥ oo, when g(n) becomes a Dirac 5-function, the transporting states are extended Bloch-like 
waves in the direction orthogonal to the field that carry the current v* pQ. 

Figure [5] depicts the results of numerical simulations for F = 0.1 and two orientations, /? = 1/3 
and (3 = {s/b — l)/4 « 0.309. (Observe that these values are half of those considered in the previous 
subsection.) As initial condition we choose the wave packet shown in the upper-left panel of Fig. [3J 
This packet is the transporting state for the field direction (r,q) = (0, 1). In his family, it is the easiest 
to construct numerically. Moreover, when F is small, it has large overlap with transporting states for 
nearby directions as well and it can be used to test the transporting regime for arbitrary field directions. 
The solid lines in the panel (a) of Fig. |6]show the wave-packet center of gravity (x(i), y(t)) for /3 = 1/3. 
It is seen that the packet moves in the direction orthogonal to the electric field: y(t) = —f3x(t), and 
that the speeds in the two directions are those predicted by semi-classical analysis. The dashed lines 
show the same quantities, now for = (VE— l)/4- The lower panel shows the wave-packets at the final 
time of numerical simulations, which appears to be still well focused. This confirms the fact that in this 
short-time, semiclassical regime, dynamics is not affected by commcnsurability of Bloch frequencies. We 
comment later on the long-time regime. 

When the initial wave-function is not in the form of a transporting state, we encounter different 
dynamical behaviors. If the initial packet overlaps significantly a transporting state, we observe a comet- 
like dynamics with the comet head moving at the drift velocity and the tail extending in the opposite 
direction. Finally, for a generic initial state with small overlap with the transporting state, the wave- 
packet dynamics is an asymmetric ballistic spreading whose dispersion increases in time approximately 
as 

a(t) « v*t/V2 . 

However, the most prominent feature of the weak field, semiclassical regime ([6]) is that it is insensitive 
to the rational versus irrational nature of the orientation /3. i.e. to the commensurability of the Bloch 
frequencies. This conclusion is consistent with the semiclassical analysis. 
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Figure 7: Dynamics of a localized wave packet for eaF/J = 2, a = 1/10 and (3 = 1/3 (solid line) and 
j3 = (\/5 — l)/4 (dashed). The lower panel shows the wave-packet dispersion in these two cases. The 
upper panel depicts x{t) and y(t) for (3 = (V5 — l)/4. 

6.2 Strong field limit 

In the large field limit, eaF/2nJ > a, the scaling law (|29l30p implies that dispersion of a wave-packet is 
inhibited for irrational field directions. Therefore, the quantum motion can only oscillate in width and 
position. On a relatively short time scale, this is also the case for "bad" rationals r/q with r, q ^S> 1, 
while ballistic spreading can be detected only for j3 = r/q with a small denominator. 

In the original frame of reference (x, y) the motion can be well described by the first momenta x(t) 
and y(t) and by the dispersion a 2 (t) = a 2 J2i m (^ 2 + m 2 )\^i,m{t)\ 2 — x 2 (t) — y 2 (t) . In Fig. [7]we display 
these data for an initial Gaussian wave packet with F = 2, a = 1/10 and for rational (3 — 1/3 and 
irrational direction (3 — (\/5 -l)/4ss 0.309 . The dashed and solid lines in the lower panel depict the 
wave-packet dispersion for these two cases, respectively. As expected, a secular increase of the dispersion 
is observed only in the rational case (3 = 1/3. In the upper panel, dashed and solid lines plot x(t) and 
y{t) in the case of the irrational field direction (3 — (\/5 — l)/4. The characteristic amplitudes and 
frequencies of oscillations of these quantities are defined by Bloch oscillations, 

X{t) = 2^ ' m = ^Fy MUJyt) ■ 

The magnetic field distorts these oscillations, the more the larger the value of a. In particular, for a 
close to its maximal value 1/2 (without any loss of generality one may consider |a| < 1/2) it becomes 
impossible to recognize Bloch oscillations in the time evolution of the first momenta. Nevertheless, the 
conclusion that a strong electric field localizes the quantum particle on a lattice remains valid. 

When derealization takes place, it is convenient to consider the rotated coordinate frame (??,£)■ This 
permits to verify numerically the scaling law (|29|30[) . In fact, we expect that the coefficient A in the 
asymptotic growth of the second moment in the r\ direction: 

M 2 (t) :=a 2 J2\Hv,0\ 2 ri 2 -^t\ 

v 

be proportional to bandwidth and hence, via Eq. ([29)) . to F 1 ~ q ~ r . In Fig.[8]we plot A versus the electric 
field amplitude F, for five cases of the ratio /3, that takes the rational values zero, one, one half, one 
third, and two thirds. Naturally, the second moment in the field direction £ is bounded in time. All data 
sets were computed as averages over different realizations of incoherent Gaussian wave-packets with the 
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Figure 8: Top: Coefficient of the asymptotic growth of the second moment in the r\ direction versus 
field amplitude F, for (r, q) = (0,1) (open squares); (r, q) = (1,1) (pluses); (r, q) = (1,2) (crosses); 
(r,q) = (1,3) (asterisks) and (r, (?) = (2,3) (full squares). Here, a = 1/10. Symbols are joined by lines 
to guide the eye. Bottom: The same, now displayed in doubly logarithmic scale. The straight lines are 
the large field estimates A ~ F 1 ^ q ^ r 
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Figure 9: (Color online) Square root of the second moment in the evolution of an incoherent wave-packet 
versus time t and field intensity F for a = 1/3 (green lines, crosses) and a = (v5 — l)/4 (red lines, 
pluses). The field is directed as (1,1). 

same initial widths. Values of A were obtained as fits over an asymptotic time range, extending to a few 
hundreds (in the units adopted) for /? = 0, and to tens of thousands, for f3 = 2/3. 

Bandwidth affects the dynamics in a second way, that might even be more relevant than the first 
in laboratory experiments. In fact, the indeterminacy principle implies that the time required for the 
dynamics to "feel" the continuous nature of the spectrum is inversely proportional to bandwidth. We 
confirmed numerically that the larger the value of r + q, the later in time the asymptotic behavior is 
achieved: the packet needs more time to unfold and to "pick up speed" : at fixed, large F, the packet 
seems to be "frozen" for times that grow exponentially in q + r — 1, see Eq. {§]). 

7 Quantum Dynamics in the general case 

In the preceding section we have considered small values of a, which insures the semiclassical analysis, 
which obviously describes the quantum motion appropriately only over a finite time scale. In this section 
we discuss the extension to the general case. 

7.1 Rational orientation, general Peierls phase 

We have already commented that, for rational orientations /3, the spectrum should be absolutely contin- 
uous for any value of the Peierls phase a. Quantum dynamics under these circumstances should reflect 
the spectral properties and be ballistic even outside the semiclassical region [13] . 

This expectation is confirmed by numerical experiments. In Fig. [9] the square root of second moment 
is shown versus time and field intensity F for two (large) values of the Peierls phase, a — 1/3 and 
a = (V5— l)/4 ~ 0.3090. The orientation is rational: /3 = 1. As usual, the incoherent packet is obtained 
averaging over realization of random phases. We observe two regions in the plot: as expected, for large 
F data are described by perturbation theory, and no difference between the rational and irrational case 
is observed. To the contrary, for small F two maxima in the ballistic speed {i.e. the slope of the linear 
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Figure 10: (Color online) Final wave-packet at time t = 5000, versus a, for v* — 0.5 (red), v* = 0.25 
(green) and v* = 0.125 (blue). Here, /? = (Vb — l)/2. Observe the "classical" sub-packet traveling in 
the positive r\ direction. 



growth of \JM-i) are observed, with a significant deep between them at around F = 1.25. Clearly, this 
feature is outside the reach of semi-classical analysis. Moreover, the smaller the value of F, the larger 
become the differences between the rational and irrational case. At null field, of course, exact analysis 
predicts ballistic motion for rational a and anomalous diffusion, i. e. quantum intermittency |14[ 115] for 
irrational a. 

To sum up, for rational j3 we find novel behaviors outside the semiclassical region only in the case of 
small electric fields. Yet, these novel behaviors are at most variations inside a general picture of ballistic 
regime. 

7.2 Irrational orientation, general Peierls phase 

Finally, we consider the case of irrational directions of the electric field. Here, the Hamiltonians (f2"U)) 
and (|2ip are not applicable, or rather they can be used in a sequence of rational approximations (r, q) to 
an irrational direction. It seems therefore that in the perturbative regime the band width (which drives 
the speed of ballistic spreading, as well as the time required to start this dynamical regime) is smaller 
than any negative power in the field intensity F. Yet, it could be a non-analytic function of this latter. 
This implies that the wave-packet is localized for large electric field at any irrational direction /?. 

For small a semiclassical analysis (insensitive to the rationality of the direction) indicates that the 
motion should be ballistic when F < F cr , i.e. F < 2iraJ/ea, a fact confirmed by numerical experiments, 
see Fig. [5] We expect this behavior to hold until deviations from the semiclassical theory emerge. In fact, 
the problem of the spectral type for any electric field intensity should be treated along the lines of [16] : 
this analysis should yield the result that the spectrum is pure point, and quantum motion is localized, 
albeit the localization length may be very large, as typically happens in two-dimensional systems [17] . 

Recall that the classical Hamiltonian (JTJ| depends only on the scaled electric field T and its orien- 
tation j3. Keeping the direction /3 fixed, as well as J — J x — J y , yields a single free classical parameter 
2 J aJ , where F is the amplitude of the electric field. It turns out that this ratio coincides with the velocity 
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Figure 11: Dispersion of the wave-packet versus time and a, for v* = 0.5. 



v* times the dimensional constant j^t (that in the units employed in this work takes the value one): 
it is therefore interesting to study the quantum dynamics at fixed v* , while varying the semiclassical 
parameter a. As a consequence, during this scan, the electric field amplitude scales as F — 2naJv*. 
Again, we observe two quite distinct dynamical regions. Firstly, for values of the classical parameter v* 
larger than one, the motion is always localized about the origin, for any non-zero value of a. 

To the contrary, a rich dynamical behavior is observed for values of v* smaller than one. In Fig. [10] 
we draw the r\ projection of the evolved wave-packet amplitude versus a, for three values of the classical 
parameter v*, and for the irrational value (3 — (Vb— l)/2. A coherent sub-packet moving in the positive 
r] direction is observed in all cases, roughly independent of the value of a, for small values of this 
parameter. Clearly, this feature should be ascribed to the classical dynamics. The region of a in which 
this sub-packet is observed diminishes when increasing v*, and disappearance of the packet is abrupt. 
For values of a larger than this critical threshold, the packet remains localized about the origin. 

In Fig. 1111 we plot the dispersion of the evolving wave-packet as a function of time and a, for v* = 0.5. 
The changing dynamical behavior for increasing a noted in the previous figures is observed also here: 
one sees an initial ballistic motion followed by a slower diffusion and by saturation. These numerical data 
provide us rough estimates of the transition time and the localization length. Notice that in this figure 
a is proportional to the field intensity F: one could therefore expect a sharp dependence of localization 
length on the external electric field, in line with the theoretical remarks presented at the beginning of 
this section. 

8 Conclusions 

In our previous work [I] we considered the wave-packet dynamics of a quantum particle in a square 2D 
lattice in the presence of a magnetic field normal to the lattice plane and an electric field, which was 
aligned along one of prime axis of the lattice (for definiteness, the y axis). In this work we have extended 
these studies to the case of arbitrary direction of the electric field vector F. We have confirmed the 
conjecture of Ref. [T] that, depending on the electric field magnitude, the system has two qualitatively 
different dynamical regimes, which we refer to as the strong and weak field regimes, respectively. 
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The new analysis has extended the validity of the theory developed in [T] and established new 
phenomena. For instance, we have used semiclassical analysis to predict ballistic derealization for small 
electric field intensities, and arbitrary field directions. Under these conditions, we have used spectral 
theory to construct non spreading wave-packets traveling at constant speed, at least for finite times. 

In fact, for weak electric field the wave-packet dynamics is governed by the cyclotron dynamics. 
This means that the packet moves in the direction orthogonal to F at the drift velocity v* , which is 
proportional to F, in close analogy with the problem of a charged particle in a free space subjected to 
crossing electric and magnetic fields. However, the presence of the lattice restricts this behavior to a 
subspace of initial conditions discussed in Sec. 16.11 For generic initial conditions we have found that the 
packet typically splits into several packets moving in the orthogonal direction with different velocities, 
both positive and negative. 

In the strong field regime the wave-packet dynamics is governed by the Bloch dynamics of a quantum 
particle in a 2D lattice. For null magnetic field, in these Bloch oscillations the packet oscillates near 
its initial position. The obvious exception to this oscillatory behavior occurs when the vector F points 
to the y direction. Here the packet spreads ballistically in the direction orthogonal to the field at a 
rate defined by the hopping matrix element. A finite magnetic field 'generalizes' this exception to the 
cases where the vector F points to a rational direction, i.e., /3 = r/q with r,q being co-prime numbers. 
However, now the rate of ballistic spreading in the direction orthogonal to the field is suppressed by a 
numerical factor proportional to (l/F)( r+9 ~ 1 ). This functional dependence implies that in practice the 
wave packet spreading can be detected only for simple rational directions with a small denominator q. 

We have found a critical field magnitude F cr , which separates in the parameter space the above 
discussed regimes, by generalizing the semiclassical approach of Ref. 1 1- This results in a strongly 
nonlinear ID classical system with quasi-periodic driving. Surprisingly, this effective system appears 
to be completely integrable in spite of the quasi-periodic character of driving, which is a rather rare 
instance from the view point of dynamical system theory. 

Finally we discussed the validity of the semiclassical approach. For irrational directions (3, at fixed 
time and classical parameter J 7 , we observe a sharp suppression of ballistic spreading, occurring when 
a overcomes a certain threshold. Combined with results for rational directions /3, this fact leads us to 
the conclusion that the energy spectrum of the system is continuous for rational (3 and pure point for 
irrational j3. This result holds for any value of the Peierls phase a, both rational and irrational. 



9 Appendix: Semiclassical Hamiltonians 



Let us obtain Eq. (fT4"| in Sec. [3] as a semiclassical approximation of Eq. (f2~Tj). We follow the same lines 
as in the derivation of Eq. (fTTj) in Sec. 12.21 As above, when \a\ is much less than one, we replace p by 
the continuous variable £ and we introduce the shift operator exp(c^). At the left hand side of Eq. (|2"T)l 
one therefore observes the action of the operator I: 



J X 

T 



exp 



lar 



N 



ide — 2nia- 



aVN 



+ h.c. 



exp 



laq 



N 



idf + 2iria- 



aVN 



h.c. 



eF(, 



Next, let us introduce the operators X — 2irat;/a and P 
[X, P] — 2ma. The previous equation becomes 



-iadc, which obey the commutation relation 



-J r cos 



P 



N 



X 



N 



Jy COS 



P 



N 



X 



N 



eaF- 
2na 



Finally, using the canonical transformation Y = (-r/VN)P + (q/VN)X, P = (q/vN)P + (r/\/N)X 
we obtain the classical Hamiltonian in Eq. (|14[) . 
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